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Abstract 

Using linear-response density-functional theory, we obtain the magnetic interactions in the several iron pnictides. The ground 
state has been found to be non-coUinear in FeSe, with a large continuum of nearly degenerate states lying very close to the 
magnetic "striped" structure. The presence of non-coUinearity also seems to be a generic feature of iron pnictides when the Fe 
moment is small. At small Rfe-Se the system is itinerant: strong frustration give rise to excess of spin entropy, long ranged 
, interactions create incommensurate orderings and strong biquadratic (or ring) couplings violate the applicability of Heisenberg 
model. There is a smooth transition to more localized behavior as Rfe-Se increases: stable magnetic orbital order develops 
which favor long range AFM stripe ordering with strongly anisotropic in-plane exchange couplings. The stabilization of the 
stripe magnetic order is accompanied by the inversion of the exchange coupling. 
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Accumulated experimental studies indicate that the recently discovered iron oxyarsenide family (Fe-As-La-0) and 
iron chalcogenide (Fe-Se) superconductors [H 0] similar itinerant antiferromagnetic metallic systems. This fact distin- 
guishes these systems from traditional high- Tc superconductors where the local-moment picture and the importance 
of strongly correlated electrons is generally accepted. While the physics of these systems can be dramatically different 
they share an important common feature: antiferromagnetic (AFM) interactions on a (nearly) 2D lattice. One can 
hope that these materials have similar physics and analysis of new metallic systems will shed some light on a nature 
of magnetic fluctuations in both type of superconductors. 

Description of magnetism in itinerant magnets is rather complicated, owing to the lack of a commonly accepted 
rigid spin picture and the disappearance of adiabaticity. The trivial use of localized models leads to oversimplification 
of the physics of itinerant magnets. As we show here, the degree of itineracy is critical to understand the ground 
states and magnetic excitations in iron pnictides. Owing to a strong competition between on-site Hund energy and 
Fe-As bonding in iron pnictides 0-11], we adopt the local density approximation (LDA) as we believe it to be the most 
appropriate readily accessible method to describe these systems [61]. 

We analyze the magnetic properties of iron pnictides, using a variety of implementations of the LSDA, including 
FLAPW, FPLMTO, LMTO-ASA. The latter includes calculations of non-coUinear configurations following the original 
prescription [8]. Full-potential results are in excellent agreement with each other, and (where they can be compared) 
very similar to those obtained in Ref. 7]. The ASA results are also in reasonably good agreement with them. For 
the linear response calculations we adopted a Green's function formalism within the ASA [3, [lol| . This technique has 
been extensively tested and applied to many different metallic and insulating systems [ol. Illlj. A similar method was 
recently applied to LaFeAsO[13. We have taken care to converge the calculations of parameters reported, e.g. using 
a 16x16x12 fc-mesh and 24 energy points for the contour integration in the complex plane. 

We begin our analysis with the low-temperature phase of the FeSeQ. At low temperature FeSe assumes an 
orthorhombic structure^, space group Cmma, with distortions of FeSe slabs similar to those observed in FeAs layers 
in the iron oxyarsenide family. Experimental lattice parameters were employed|2i] with Fe at (4b) and As at (4g) 
Wyckoff positions. In both FeSe and LaFeAsO, there is an internal parameter z not fixed by symmetry, corresponding 
to the height of the chalcogen above the (possibly distorted) square of Fe atoms lying in the xy plane, z fixes the 
Fe-Se distance RFe-Se, but we treat it here as a parameter because of the ambiguity in determining it even while the 
magnetic properties depend strongly on Rf^e-Se(3|. The similar dependencies of magnetic moment on volume have 
been observed in other iron based magnets (see Ref. Ts*] for the review) . 

Our initial magnetic structure was obtained by relaxing a 32 Fe atom supercell (with z=0.25). Spins were initially 
put in random (noncoUinear) orientations and their orientations relaxed along with the spin and charge densities, 
as described in Ref. 14]. In this way several magnetic structures have been identified (Fig.l). The self-consistent 
(minimum-energy) coUinear solution was found to be an AFM configuration in the striped geometry, with a Fe 
moment M ~ 0.85/iB. As Fig. 1 shows, M is very sensitive to Kpe-Se- The Neel state is the least stable among the 
structures indicated in Fig. 1: it is '^l.O mRy/Fe less stable than the striped structure (and nonmagnetic at the LSDA 
equilibrium geometry). All structures exhibit similar dependence on M{Rpe-Se)- At very small moments {M<0.2ij,b) 
orderings {b),{c),{d) in Fig. 1 have lower energy then the striped configuration (a). For M>0.95/iB (large Kpe-Se) 
the dimer structure (Fig.(le)) has lowest energy. Such sequence of ground states clearly indicate the significance of 
couplings beyond first two nearest neighbors. 

The following picture emerges from studies of the pairwise exchange coupling parameters Jij and non-magnetic 
susceptibility Xij hi different structures, as a function ofUpe-Se- We denote xoo as the on-site susceptibility and xoi 
and xo2 that of the two nearest neighbors (NN). Much of the behavior we observe depends on these three quantities. 
Define 77 = X02/X01 as the "magnetic frustration" parameter (see below). For small M (small Kpe-Se), Xoo is rather 
small ('^14-15 Ry~^). Noting that the Stoner parameter I is JwO.068 Ry for Fe, the local moment criterion Ixoo > 1 is 
marginally, or not quite satisfied (see Ref.[l0[). Thus the system already at this stage can be classified as a marginally 
itinerant system: that is the magnetic moment depends on the environment for its stability as opposed to localized 
behavior, where the local moment is large and weakly dependent on environment. We find both xoi and X02 are large 
and negative. So, the criterion for AFM pair formation I{xoo — Xoi) > 1 along (1,0), (0,1) and (1,1) is fulfilled while 
the criterion of FM pair formation I{xoo + Xoi) < 1 is not. FM fluctuations are thus strongly suppressed while AFM 
interactions arc frustrated, because both xoi < and xo2 < 0, with rj ~ 0.4. Increasing Kpe-Se (or volume) from a 
small value, a magnetic moment appears as a consequence of the sum of exchange fields from all sites, supporting the 
itinerant nature of the ground state there. This is characterized by a = Jx((7=0) /xoo ~ 1] 1 which differentiates localized 
and itinerant svstems[18j. For small M, a^O.S. As Ri?e-Se increases to the magnetic instability point Ix — ^, Ix (q) 
approaches 1 in the small region of wave vectors corresponding to striped and dimer orderings. At low moments AFM 
states can be created and destroyed by arbitrarily weak interactions. 

As Kpe-Se increases further, M increases smoothly and 77 decreases, reaching ~0.2 for large M. Thus the magnetic 
state evolves to a less itinerant and less frustrated condition and even FM fluctuations can exist, a decreases smoothly 
with M, but the magnetic behavior continues to be borderline between itinerant and localized for intermediate R_Fe_se, 
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and the frustrated nature of a weakly magnetic state is conserved for all reasonable distances 'R.pe-Se- Only in the 
strongly magnetic state is the situation drastically changed. 

In Table 1 we present both longitudinal and transverse components of the pairwise magnetic "exchange" parameters 
Jij for the striped structure with M ~ O.Qfj^B (near the theoretically optimized LSDA value). Both components of 
the exchange matrix J have similar amplitudes, indicating that both longitudinal and transverse spin fluctuations 
must be taken into account. The NN coupling between atoms with parallel spins (along x, r=l/2,0,0) is negative 
("frustrated") while the exchange between pairs with antiparallel spins (along y, t=(0, 1/2,0)) is positive (also for 
T=(l/2, 1/2,0)). For more distant (1-5 to 1-9) in-plane neighbors all parameters have similar amplitudes with different 
signs. While two of the NN parameters indicate a stable pair configuration. Fig. 2 and the Table show the presence 
of rather unstable FM pair alignments in plane, and weak interaction along the z-axis (see Jqoi in the Table 1), 
making this system very nearly a 2D AFM. However, even such weak interplanar coupling may lead to 3D long-range 
order and finite temperature transition. The insert of Fig. 3 shows how the Fe moment in the striped structure 
depends on anisotropy of the exchange coupling of NN atoms. The origin of this anisotropy is as follows. For the 
Neel structure both Ji^ = Jioo and Jiy = Jqio are positive and equal (reflecting symmetry) while J2 = Jno is 
negative, reflecting instability of FM alignment along (1,1). For the striped structure at small R_Fe-Se ("bad" nesting 
conditions) the system still reflects the frustrated Ji- J2-J3 model-like behavior and our Jix and Jiy have opposite 
signs with similar amplitudes (Fig. 2). As M increases Ji becomes anisotropic, reflecting a broken symmetry induced 
by magnetic analog of Jahn- Teller effect [TH. [l6l| . For M ~ 1.5^b the anisotropy is very large and soon 'frustrated' 
exchange coupling along FM line disappears and the effective sign of this interaction becomes positive. Thus, large 
moment state does correspond to localized Heisenberg model and is absolutely stable with no sign of antiferromagnetic 
frustrations. At high temperature both x and y directions are equivalent, but as T decreases the barrier between x 
and y orientations will prevent fluctuations between them and the spin symmetry-breaking phase — with attendant 
structural distortions-will be frozen in. Fig. (Ic) depicts the "rotator" structure — a state intermediate between x- 
and y- aligned striped phases. At M~ I^b the LSDA energy of this state is 0.8 mRy/Fe relative to the striped state. 

The spin wave velocities are not formed merely by NN interactions in such a system, and the contribution from 
distant Fe atoms (3-4 A) with antiparallel spins is signiflcant, as expected. We found significant in-plane and out- 
of-plane anisotropics of spin velocities. For instance Vx/'^y = 1-6 for M — l.OlfiB- This anisotropy has not been 
determined experimentally. The calculated energy gap in the spin wave spectrum due to spin orbital coupling [l3| 
is A = 0.15 mRy. At larger q we see deviations from the adiabatic magnon picture. However the adiabaticity 
parameter [l^ for I^b is ~ 0.2; so we believe the system is already less itinerant and more adiabatic, therefore 
antiferromagnons can be defined not only at small q. However longitudinal fluctuations of similar strength are always 
present and should be seen in experiment. These features of spin excitation spectra will be reflected in the temperature 
dependence of the sublattice magnetization and will modify the enhancement of the linear speciflc heat. Also at those 
Upe-Se where the spin fluctuations are very localized around a particular Q one may expect weak temperature 
dependence of the susceptibility. 

Because parameters Jij strongly depend on magnetic configuration (which may not be a ground state) and on 
symmetry-breaking, it is crucial to analyze the magnetic interactions for several magnetic orderings. The universal 
picture that emerges is that first and second NN interactions (both transverse and longitudinal) on the square Fe 
sublattice are AFM and very strong, with their absolute values and the ratio ry depending on Kp^-se- The role of 
longitudinal fiuctuations is crucial at larger M . Even so the relatively small more distant neighbor couplings (and 
biquadratic terms) are not negligible and important for the appearance of non collinear fiuctuations at small moments 
as we now describe. 

To check the stability of our collinear " ground" striped state and avoid uncertainties connected with the long- wave 
approximation [1^, we performed spin dynamics simulations with small deviations between FM coupled atoms keeping 
first or second NN ordering antiferromagnetic. These excitations include ji and j2 modes (Fig.l) and fictitious spin 
spiral (SS) orderings[8]. Our results demonstrate that for those Rpe-Se with M< I^b (Fig. 3) spirals (Q -I- (0,g,0)) 
which rotate AFM planes against each other exist only in some limited spin phase space around Q of the striped 
structure, clearly establishing the existence of the magnetic short range order (MSRO). This MSRO is significant and 
can exist far above the critical temperature. However, the phase volume of such long-wave excitations is rather small. 
As Kpe-Se increases the short ranged excitations with much larger volume of phase space become available and it 
leads to stabilization of the different magnetic ground state or changes the character of excitations. 

The total energy of SS configurations which rotate FM spins along the stripes, i.e. (Q + (^,0,0)) appears very 
different. Fig. 3 shows the energy is nearly independent of q for a significant region of q. Such a flat surface in the 
itinerant AFM usually is a precursor of frustrated interactions. This is exactly the present case, as evident from 
the negative value of Jix (Fig 2). Moreover SS calculations show that the collinear striped structure {q=0), which 
has strong AFM coupling between stripes and no staggered magnetization at flnite q, is not the minimum-energy 
state, although the stabilization of this particular SS is very small (~ 0.01 mRy) with SS q consistent with the 
amplitude of J3, J4 couplings. We also found that ji fluctuating mode (which has staggered magnetization) has very 
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low excitation energy, so states in this fluctuating direction are practically degenerate and only starting at 15-20° do 
we see a small energy increase ^ 0.03 niRy. At larger angles the energy rapidly increases, so the "zig-zag" (45°) state 
has energy ~ 0.3 niRy/Fe. The nematic j2 mode has much higher energy at large moments (which quite natural for 
the biquadratic or ring type of exchange coupling), while at small moments it has lower energy than the ji mode, 
supporting the increase of frustration parameter 77. Both "zig-zag" and "rotator" structures are indicators of the 
strength of biquadratic exchange (K) and its importance for the phase formation and spin wave spectrum gap. It is 
unusual to have biquadratic exchange as large as half of bilinear exchange parameter for the system with relatively 
small values of effective spin value (1/2 or 1). However, we obtained a similar behaviour in many other pnictides 
and believe that this is a generic feature of these systems. This fact, in turn, may create favorable conditions for the 
applicability of model, suggested in Ref. [l3| • The biggest difference with that model consists in the origin of barrier 
energy. In our calculations this barrier exists already on adiabatic energy surface as a result of complicated spin-spin 
interactions in the ground state, providing even stronger support of conclusions of Ref. . In addition, the low 
energy of ji and corresponding 'nematic' type of spin ordering (j2) can be responsible for the disappearance of LRO 
in plane due to developing at finite temperatures fluctuations between x and y stripes. This mechanism of magnetic 
transition can be considered as an alternative to usual formation of two-dimensional state which is widely used in 
high temperature superconductors study. Besides being important for the magnetic transition, these fluctuations must 
appear at the finite temperatures, may persist in the paramagnetic case and should be observable in the experiment. 
These type of excitations are very unusual: they do not exist at low temperatures in the ordered state, but will appear 
closer to the Neel temperature or above it. 

Thus, our exchange coupling and SS calculations indicate the presence of two main competing mechanisms: one 
is an instability (driving force) towards formation of AFM order ("good" nesting conditions), and second frustration 
related to the strong negative couplings of first and second NN on a square lattice (with rj ~ 0.2 — 0.4). While in 
the J1-J2 model Ts'l rj is the only frustration parameter, our calculations reflect a more complex picture: the relative 
strength of all frustrations (as any interatomic coupling) is also determined by the "itineracy" characterized by a. 

The latter mechanism gives rise to large spin fluctuations related to J3 or J4 pair interactions and is particularly 
important at small M (itinerant limit) with strong MSRO and limited available phase space with long-wavelength 
excitations (analogous non-collinearity we also found in CaFe9As9^1(l|). These fluctuations can even destroy the long 
range AFM order at lower Kpe-Se (or lower volume or c/a ratio). For larger Kpe-Se improving nesting conditions 
and a naturally developing unstable longitudinal mode enhance the first effect: a stable, ordered AFM state forms 
with large more localized M. The itineracy parameter a coincides with the "adiabaticity" parameter of Ref.fTsj 
and indicates that strongly frustrated metallic systems simultaneously have strong non-adiabatic effects. Thus, our 
calculations shown that transverse frustrations are rapidly decreasing as the system becomes more localized. 

To illustrate our result for FeSe we map the results obtained above onto Heisenberg Ji- J2-J3 modeljlBl. Fig. 4 clearly 
shows how relatively weak exchange coupling between third NN transforms the collinear ground state. Especially 
important are the third and forth NN coupling when Jiw2 J2, when the non-collinear state of different symmetry can 
be stabilized by a very small coupling beyond second NN. In strongly magnetic case (M=l/is) when NN coupling 
appear to be very anisotropic the ratio between main parameters of spin Hamiltonian are following (122 family): Jij, 
« 30 — iOmeV, Jix ~ 0.2Jij,, J2 « 0.7Jij^, Jz ~ 0.15Jij^, Ki w (0.2 — 0.3)Jia and magnetic anisotropy ^0.1 nieV. 

Despite theory predicting very similar properties for many families of iron pnictides, experimental evidence for long 
range magnetic order in the normal state have not been found in all compounds. In particular, no long range order 
of any kind been found in LiFeAs or FeSe, nor have local moments been found so far. In Tab. 2 we summarize the 
calculated exchange coupling obtained for the different iron pnictide systems, including some with no experimentally 
established LRO. All these systems clearly demonstrate a stable AFM stripe order at large moments. The exchange 
coupling shows rather itinerant behavior for smaller moments (large range of parameters), and more localized for 
the larger moment case. Due to large anisotropy between first NN exchange parameters the striped structure is 
very stable and the isotropic J1-J2 Heisenberg model is totally not applicable for iron pnictide systems at larger 
magnetic moments where this anisotropy is very large. However, this model might be relevant at smaller moments 
when Ji is close to J2 . The similarity in the exchange coupling and the comparable stability of the stripe order clearly 
contradicts the available experimental data for these systems. Assuming that more experiments with better single 
crystals confirm that LRO in LiFeAs and FeSe is absent, the fact that theory contradicts these observations must be 
related to LSDA shortcomings. In particular, magnetic zero-point fiuctuations (absent in LSDA) may well be strong 
enough to suppress the moment. Alternatively, there may be problems even in the LDA description of bonding in 
iron pnictides at the independent-particle level. 

Finally, we calculate the imaginary part of the bare dynamical transverse spin susceptibility (q? ^) i which 

provides insight into Stoner excitations [l^. Calculations were performed for the Neel structure at 2=0.25, which 
corresponds to a magnetic moment M=1.2fj,B- Along (110), x"^ (Qi"^) is particularly interesting; see Fig. 5a. 
Stoner excitations appear even for uj—^-O for small q. As q increases they disappear completely; see q = (0.25,0.25). 
Finally, they reappear at still higher q. 
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Most of these novel features can be explained in terms of the unique nature of the Fermi surface. The Fermi surfaces 
of FeSe (and the Fe pnictides) are approximately cylindrical, and are centered either at F or at the (tt, tt) point in the 
xy plane, as shown in Fig. 3 of Ref.[7] Cylinders at F are of hole character; those at the zone corner are of electron 
character, q is the connecting vector that links an electron at k + q to a hole at k. x°~' (<li'^) is obtained from 
a sum over all allowed k (i.e. any k for which an electron state at k + q and a hole state at k are separated by an 
energy difference uj). We have identified the following facts, which can explain the behavior in Fig. 5a. 

(a) For q < (0.25, 0.25, 0), k points only couple to points k + q in the same pocket. 

(6) For q > (0.25, 0.25, 0), k points only couple to points k + q in other pocket. 

(c) For q = (0.25, 0.25, 0), no two k points can be connected. 

(d) For q > (0.375,0.375,0). there is a strong inter-pocket excitation. 

Along (100) the story is different. Since the two pockets are separated by (1/2, 1/2, 0), they can not communicate 
for q=((700) (at least at low energy). So only intra-pocket excitations, F— F, M— M, are allowed. This picture is borne 
out in Imx°"' (q, w) , shown in Fig. 5b. Imx°~' (Qjw) is already present for small w, e.g. 50 meV, when q is small; 
but almost disappears when q reaches 0.25, and completely disappears for larger q. 

While we have not yet attempted to calculate spin fluctuations, the fact that large Stoner excitations are present in 
a significant fraction of the Brillouin zone offer strong indications that low energy spin-flip excitations are important. 
They can strongly affect the single-particle description of the magnetism; in particular they reduce the magnetic 
moment M . These fiuctuations may explain the discrepancy between observed magnetic order, and that predicted by 
the LDA in case of FeSe. 

In conclusion, two competing mechanisms in the newly discovered superconductors have been identified: strong 
geometric frustrations dominant for small moments, which allow distant interactions lead to a symmetry breaking, 
forming a continuum of non-collinear states. As M increases, magnetic Jahn- Teller effect [l^ creates better nesting 
conditions with stabilzation of an AFM striped configuration. We have found third and fourth neighbor couplings J3 
and J4 to be the driving force for non-collinear spin fluctuations. We believe that is the generic feature of iron pnictides 
when the Fe moment is small. The calculated exchange coupling and the corresponding spin wave excitations along 
z-direction are relatively small, can be easily destroyed by temperature and two-dimensional character of magnetism is 
likely to be realized at higher temperatures or paramagnetic state. We also find that strong frustration — a condition 
necessary for excess entropy in order to admit a superconducting state — are favored in the itinerant magnets when 
conditions for adiabaticity and long wavelength approximation are not well satisfied. Here longitudinal fluctuations 
make an additional contribution to the spin entropy. As moment increases, the in-plane exchange coupling along 
'ferromagnetic' line changes its sign, so the stripe order becomes absolutely stable for such 'no frustration' localized 
magnetic moment case. We identified those magnetic excitations in pnictides that may exist below Neel temperature 
and contribute to observable properties. Together with traditional spin waves, in this system coexist Stoner type 
of excitations (electron-hole pairs) and longitudinal spin fiuctuations, especially at small moment. The presence of 
low energy Stoner excitations ultimately identifies this system as itinerant magnetic system. In addition to such 
traditional excitations, this system has very peculiar four fold symmetry planar spin fluctuations with the energy of 
15 meV and below (modes (b)-(d) on Fig.l) with some of them being similar to 'nematic' spin order. Our results 
indicate that the localized spin description is somewhat academic for this system if magnetic moment does not exceed 
1/is. Overall, the magnetic transformation in iron pnictides can be described as a transition from frustrated itinerant 
moment phase to symmetry broken orbitally ordered localized magnetic moment phase with "striped" order. The 
transition happened around M—lfis- 

Since this paper has been submitted new neutron experimental data have been published for FeTei_a;Sea;^3|- The 
authors observed the presence of two-dimensional incommensurate excitations which may reflect the non-collinear 
ordering identifled above. 
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Table. 1 The first eight Fe-Fe transverse Jt, and longitudinal Jl exchange parameters in FeSe. Column n denotes 
the number of equivalent nearest neighbors, together with the spin orientation. R denotes the distance to the central 
atom and 7" the connecting vector in units of the lattice parameter a. Jt is related to the Hcisenbcrg model J as 
JSjSj = Jt and can be defined as a second derivative of the LSDA total energy with respect to moment variation. 
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CaFe2As2 


1.61 


-1.5 
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-12.1 
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1.07 
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-13.5 


0.8 


-0.1 


4.5 


-2.2 


-0.3 
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Table. 2 The several Fe-Fe Heisenberg model parameters J in different families of iron pnictides (in meV). These 
parameters have been obtained from theoretically calculated Jt by normahzing it by scalar product SiSj where 
Si=M,/2. 

Fig.l The dependence of magnetic moment on Fe site as a function Fe-Se distance (Rpe-Se) for the striped (lowest 
energy) and simple Neel (highest energy) magnetic structures. Also shown are schematics of the magnetic structures 
considered in the calculations: (a) striped structure (b) "zig-zag" structure. A transition from a striped to "zig- 
zag" structure (a 'ncmatic' type of state) leaves AFM ordering between first NN atoms unchanged (ji mode), (c) 
"rotator" structure (or so called 2k structure). Here AFM ordering between second NN is unchanged (j2 mode) (d) 
mixed structure with AFM moments along diagonal and non magnetic NN atoms (e) dimer structure (f) simple Neel 
structure. 

Fig. 2 Magnetic interaction parameters in FeSe for the striped structure, as a function of Fe magnetic moment M 
(controlled by varying Kpe-Se)- Jix = Jioo and Jiy = Jqio are the NN interaction along x and y respectively; 
J2 = Jiio is the second-NN interaction in the plane (vectors are shown in the Table.) Interactions between planes 
are small (~ 6% of J2). Insert shows the anisotropy of the NN exchange coupling as a function of Fe moment. The 
inversion of the sign of the exchange coupling paramer Ji^ occurs at Mssl.88 ^b- 

Fig. 3 The total energy of a spin spiral as a function of spiral wave vector q near the AFM striped structure. Top: 
SS for q=(0,q,0)27r/a. Near ^=0.25 (not shown), the moment disappears entirely. Bottom: SS along the stripe, for 
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q=(g,0,0)27r/a. The minimum energy occurs at 5 w 0.02, corresponding to a slow spiral along "frustrated" direction 

X. 

Fig. 4 The phase diagram of J1-J2-J3 model and calculated non-magnetic susceptibilities ratio for FeSe (central red 
star). Other stars correspond to 2% and 4% of hole (right) and electron (left) doping of FeSe. In all cases the finite 
value of J3 coupling leads to the stabilization of non-coUinear order of different symmetry. The relative strength of 
J3 is increased as moment decreases and the system becomes more itinerant. 

Fig. 5 Bare transverse spin susceptibility (q, w) as a function of co for different q points along (100) and (110). 
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